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ABSTRACT 

The evolution of the phase-space density profile in dark matter (DM) halos is investigated by 
means of constrained simulations, designed to control the merging history of a given DM halo. Halos 
evolve through a series of quiescent phases of a slow accretion intermitted by violent events of major 
mergers. In the quiescent phases the density of the halo closely follows the NFW profile and the 
phase-space density profile, Q(r), is given by the Taylor & Navarro power law, r _/3 , where (5 w 1.9 
and stays remarkably stable over the Hubble time. Expressing the phase-space density by the NFW 
parameters, Q(r) = Q s {r / R s )~@ , the evolution of Q is determined by Q s . We have found that the 
effective mass surface density within i? s , E s = p s R s , remains constant throughout the evolution of a 

— 5/2 

given DM halo along the main branch of its merging tree. This invariance entails that Q s oc R s and 

^ i*2 5/2 ( \ ^ 

Q(r) oc S s ' R s ' [r/R s \ . It follows that the phase-space density remains constant, in the sense 

— 5/2 

of Q s = const., in the quiescent phases and it decreases as R s in the violent ones. The physical 
origin of the NFW density profile and the phase-space density power law is still unknown. Yet, the 
numerical experiments show that halos recover these relations after the violent phases. The major 

5 /2 

mergers drive R s to increase and Q s to decrease discontinuously while keeping Q s x R s = const. 
The virial equilibrium in the quiescent phases implies that a DM halos evolves along a sequence of 
NFW profiles with constant energy per unit volume (i.e., pressure) within R s . 

Subject headings: cosmology: dark matter — galaxies: evolution — galaxies: formation — galaxies: 
halos — galaxies: interactions — galaxies: kinematics and dynamics 
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1. INTRODUCTION 

The dynamics of dark matter (DM) halos in the Cold 
Dark Matter (CDM) cosmology can be easily formulated 
as the classical Newtonian TV-body problem. Yet, the un- 
derstanding of the equilibrium configuration of the DM 
halos defies a simple analytical approach. The lack of 
analytical understanding is often compensated for by 
numerical simulations which provide empirical knowl- 
edge. The cumulative work in cosmology over the last 
decade or so has led to a broad consensus about two ba- 
sic facts that describe the equilibrium structure of DM 
halos. One is that the spherically-averaged density pro- 
file p(r) is well approxim a ted by the so-called NFW pro- 
file ([Navarro et all 119961. fl997|) or some close variants 
of it ([Moore et alJl!998t Uing fe SutJl2000l : IKlvpin et all 
120011) . The other is the power law behavior of the phase- 
space density profile, namely Q(r) — p(r)/er 3 (r) oc r^ 13 , 
with (3 ~ 1.9, where a(r) is the velocity dispersion 
(|Tavlor fc Navarrdl200ll ). 

Two seemingly orthogonal approaches to the prob- 
lem of the origin of the equilibrium structure in 
DM halos exist. One assumes a monolithic col- 
lapse of a halo that can be approxima ted by the 
spherical infall model (jGunn fc Gott) 119721 ) . The ap- 
plication of the model to the cosmological context, 
where the shell crossing has to be explicitly ac- 
counted for, has resulted in the so-called secondary 
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Rvden fc Gunnlll987t IzTaroubi fc Hoffman! 1 19931: iNusser. 
20011: lLokas fc Hoffmanll2000D . The SIM has been tested 



against the TV-body simulations and has proven to faith- 
fully reproduce the density profile of simula t ed DM halos 
(iQuinn et all 119861: lEfstathiou et al J 1 19881 ICrone et all 
119941: lAscasibar et alj|2004t l2007t ). A closely related vari- 
ant of the SIM replaces its dependence on the primor- 
dial over-density of the proto-ha lo by the mass accre- 
tion history (MAH) of the halo (INusser fc Shethl 119991: 
ILu et alJl200a ISalvador-Sole et alJl2007l) . The SIM and 
its MAH variant can reproduce also the power law be- 
havior of the phase-space d ensity (|Austin et al.l 120051 : 
iGonzalez-Casado et"aT1l2007t ). 

However, a close inspection of the iV-body simu- 
lations reveals that a DM halo evolves very differ- 
ently from a monolithic quasi-spherical collapse. In 
fact halos are numerically observed to go through a 
sequence of mergers, some labeled as major merg- 
ers in which the two main progenitors are of a sim- 
ilar mass, le ading to eme r gence of the NFW den- 
sity profile (ISver fc Wh itel jl99ll IDekeletal l [2003 



Subramanian et al. 1 120001) . |R omano-Diaz et al.l (|2006l : 
20071 hereafter Paper I and II) studied the formation 



and equilibrium configuration of halos by means of con- 
trolled TV-body simulations, with the initial conditions 
set by constrained realizations of Gaussian fields. These 
simulations were designed to address issues of how the 
merging history affects the DM halos. The emerging pic- 
ture is that of a halo evolving via a sequence of quiescent 
phases of a slow mass accretion intermitted by violent 
episodes of major mergers. In the quiescent phases, the 
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density is well fitted by an NFW profile, the inner (NFW) 
scale radius R s and the mass enclosed within it (M s ) re- 
main constant, and the virial radius (i? V ir) grows linearly 
with the expansion parameter (a). In the violent phases, 
the halos are not in a dynamical equilibrium, but are 
rather in a transition state, resulting in a discontinuous 
growth of R s and i? v i r . In such a picture a halo is defined 
in the context of a merger tree — at any given time it is 
taken as the most massive progenitor along the branch 
leading to the final halo. 

A direct comparison between the SIM and nu- 
merical simulat i ons h as been conducted recently by 
lAscasibar et al.l (|2007f) . This comparison is based on 
selecting the DM halos from a cosmological simulation, 
tracing them back in time, and recovering their individ- 
ual initial conditions. The SIM has been applied to the 
'primordial' density profiles, and their virial density pro- 
files have been calculated for different redshifts. The SIM 
calculated profiles provided a good match to the evolu- 
tion and structure of the simulated clusters. This is en- 
couraging. However, a rigorous fundamental theory that 
can accommodate both the spherical monolithic collapse 
and the major merger-driven evolution exhibited by the 
simulations is still missing. This motivates us to look 
further into the phenomenology of the phase-space den- 
sity, to gain a further insight into this seemingly simple, 
yet complicated, problem. 

The study of the DM halo evolution has been heav- 
ily focused on the density profile, while the evolution 
of the phase-space density has b een largely ignored. 
iPeirani fc de Freitas Pachecol (|2007f ) presented one of the 
few studies of Q(r) evolution. They found that Q, de- 
fined as a global quantity characterizing the halo as a 
whole, is generally decreasing with time. Specifically, it 
exhibits an early rapid decrease (at redshifts z > 6.5) and 
a late slow decrease. Here we aim at studying the evo- 
lution of the phase-space density within the framework 
of the NFW scaling and the dynamical picture formu- 
lated in Papers I and II. In particular, we shall rely on 
the empirical fact that the dynamics of the DM halos is 
constrained by a new invariant of motion. This invari- 
ant tags a halo along the main branch of the merging 
tree. The NFW scaling and this invariant of motion pro- 
vide a full description of the cosmological evolution of 
the phase-space density for individual halos. 

The structure of the paper is as follows. The analysis 
of the DM halos in Papers I and II is briefly reviewed in 
General considerations of the evolution of DM halos 
arc given infj3l The general evolution of the DM phase- 
space density profile is described in 2] and self-similarity 
and scaling relations are given in [J5j A general discussion 
follows in SJ5] 

2. NUMERICAL EXPERIMENTS 

In Papers I and II, we investigated the cosmological 
evolution and structure of five DM halos by setting the 
initial conditions of the simulations using constrained re- 
alizations of Gaussian fields. Our basic motivation was 
to perform controlled numerical experiments designed to 
study the dependence of the evolution and structure of a 
iven halo on its merging history. The lHoffman fc Ribakl 
199 If ) algorithm of constrained realizations of Gaussian 
fields has been used to set up the initial conditions. Pa- 
pers I and II present the analysis of five different mod- 



els of a given DM halo, evolving along various merging 
histories. The models were simulated within the frame- 
work of the open CDM (OCDM) model with Q, = 0.3, 
h = 0.7 and as = 0.9, where fig is the current cosmologi- 
cal matter density parameter and a% is the variance of the 
density field convolved with a top-hat window of radius 
8/i _1 Mpc used to normalize the power spectrum. This 
model is very close to the 'concordance' ACDM model 
in its dynamical properties. The models are labeled as 
OCDMa, OCDMb, etc., following the notations of Pa- 
pers I and II. Here we add a new model, run within 
the flat-A cosmology wi th the parameters of the WMAP 
three years data base dSpergel et al.l I2006T ). The con- 
straints used to set the model are similar to OCDMa 
and the model is labeled as WMAP3a, yet it is performed 
within a 256 3 computational box and it starts from a dif- 
ferent random realization of the initial conditions. A full 
description of the numerical simulations, the application 
of the algorithm of constrained realizations and the nu- 
merical code are given in Papers I and II. The WMAP3a 
model is one in a series of runs in the WMAP3 cosmology, 
to be reported in a forthcoming paper. The NFW fitting 
algorithm is described in Paper II and the phase-space 
density power law profile is fitted in a similar way. 

3. GENERAL CONSIDERATIONS 

The density profile of DM halos is well approximated 
by the NFW profile, 

4p s 



p(r) 



(r/R s )(l + r/Rs)^ 



(1) 



in which the characteristic density (p s ) and scale radius 
(R s ) define the NFW profile. Defining the halo as a col- 
lection of particles in a spherical (say) volume in which 
the mean density equals some critical over-density (which 
is in general redshift dependent) times the mean cosmo- 
logical density, the virial mass (M vir ) and radius (R V i r ) 
of the halo are determined. It follows that i? v ir and R s 
(or equivalcntly M v - lr and M s ) are the two independent 
parameters that define an NFW halo. 

Assuming the NFW parameterization, we can write the 
phase-space density profile as 



Q(r) = Q s Q 



R, 



where 



Q(x) 



P{x)/Ps 

[a(x)/a s 



(2) 



(3) 



(4) 



a s is the mean velocity dispersion within R s and x = 
r/R s . To the extent that the DM halos are fitted by the 
NFW profile, their Q(x) profile should obey a universal 
relation. Their cosmological evolution is then determined 
by the evolution of Q s . 

4. EVOLUTION OF DARK MATTER HALOS 

The cosmological evolution of the main halo of the six 
different models is best presented by Fig. [IJ which shows 
the evolution of R s and Q s . The halo goes through vi- 
olent episodes of major mergers and quiescent phases of 
slow accretion. In the OCDM models the violent events 
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are well separated by the quiescent phases characterized 
by an NFW structure (Papers I and II). The WMAP3 
model goes through an early phase of successive violent 
mergers, frequent enough so that the halo does not re- 
lax to an NFW-like configuration in between. This early 
phase is followed by a quiescent phase, which is slightly 
perturbed by mergers not strong enough to be qualified 
as major. This is clearly shown by the behavior R s and 
Q s , where R s (Q s ) increases (decreases) discontinuously 
in the violent phases and remains constant in the quies- 
cent ones. One should note that in the violent episodes 
the halos are not in an equilibrium and therefore the 
NFW fitting is very unstable and the resulting R s and 
Q s parameters are quite erratic. This is reflected by the 
spiky behavior of these quantities in the violent phases. 

Evolution of the dimcnsionless phase-space density 
profile, Q(x), is presented in Fig. [21 where the Q(r/R s ) 
profiles of the OCDMc and WMAP3 models are evalu- 
ated at different epochs, covering the time interval from 
z = 5.3 (z = 3.35) for the OCDMc (WMA3) model to 
the present epoch. The evolution of the other models is 
virtually identical to the ones shown here. The Q{x) pro- 
files are very closely approximated by a power law, with 
a fractional deviation of less than a twenty percent (lower 
panel of Fig. [2]). Fig. [2] shows not only the power law 
nature of the profiles but also that indeed the Q s scaling 
renders the Q(x) profile to a universal time independent 
power law. The evolution of the exponents of the Q(x) 
fitted power law of all the models is shown in Fig. [3] (3 
displays a bumpy behavior, with bumps corresponding 
to the major mergers, but generally staying within the 
ps 1.9 ± 0.1 range. Overall, given the violent character 
of the halo's evolution, the robustness of the (3 range is 
remarkable. 

A close inspection of the halo evolution reveals that 
the product S s = p s R s remains approximately invariant 
as the halo evolves along the main branch of its merging 
tree (Fig. 0J. The value of S s fluctuates around its mean 
value in two different modes., undergoes fast and corre- 
lated fluctuations of a small amplitude (see next section 
and Paper II). In the violent episodes it exhibits large 
deviations from the mean value. This behavior is asso- 
ciated with timings of major mergers when halos are far 
from the equilibrium. The NFW fitting fails here and 
the resulting NFW parameters are ill-defined. Ignoring 
the spikes of the violent phases, and averaging over the 
jittery fluctuations in the quiescent phases, E s retains its 
value along the halo evolution. While the invariance of 
E s in a given quiescent phase is not surprising, its ability 
to retain the value in and across the major merger events 
is not obviously expected. 

5. SELF-SIMILARITY AND SCALING RELATIONS 

Assuming the empirical finding of S s invariance, we 
shall study its ramifications for the evolution of the 
phase-space density. In the quiescent phases, in which 
the density follows the NFW profile, a halo is expected 
to be in virial equilibrium, as corrobrated by Fig. 13 of 
Paper II. In particular the virial ratio, evaluated within 
R s , should have a constant value, namely, 

f — ps const. (5) 



The virial ratio differs from unity because the inner part 
of the halo does not constitute an isolated system, and 
its value depends on the shape of the density profile. 

Assuming the constancy of the virial ratio within R s 
and the invariance of S s , 

a^oc^ioc^ = S 2 (6) 

The evolution of the phase-space density profile in a given 
halo is described by 

QWocS- 1 /^-^^)^. (7 ) 

Consequently, for a given halo the following invariance 
holds: 

Q S R 5 S /2 PS const. (8) 

This prediction has been tested against the five models 
of Papers I and II and the ACDM halo. Fig. [5] shows the 

cosmological evolution of Q s R b J 2 of these models. As 
with all other quantities that characterize the DM halos, 

5/2 

the product Q s Rs shows a jittery behavior in the quies- 
cent phases and strong fluctuations in the violent phases, 
in which the NFW parameters are ill defined. Apart from 
this, it remains constant throughout the evolution. Only 
Model C exhibits a small deviation from this invariance. 
In its first quiescent phase Q s rV 2 is larger by a factor 
of < 1.5 than its asymptotic value. 

The p s R s invariance is to be distinguished from the 
p s oc Rj m scaling found in its t he low-ampl i tude j itter 
(Paper II) and from M s oc R% of IZhao et all (|2003t ). In 
Paper II we show that m ps 1.39 in the last quiescent 
phases of various halos, while at earl y times m p s 1.59. 
This stands in good agreement with IZhao et aLl 's a ps 
1.44 in the "slow accretion phase" and a ps 1.92 in the 
"rapid accretion phase." The m ps 1.39 scaling is ob- 
tained by analyzing each quiescent phase, and it reflects 
the fluctuations of p s and R s around their mean values 
within that phase. This correlation appear to be driven 
by density fluctuations that originate in the region be- 
tween the cusp and R s , where the density slope varies 
between -1 and -2. The mean values of p s and R s change 
from one quiescent phase to another. The m ps 1.59 
reflects the variation across the different phases. The as- 
sociation of the m ps 1.59 and a ps 1.92 with the early 
times of the evolution of halos stems from the fact that 
the violent phases are more abundant at early times. The 
combined analysis of Paper II and the present work yields 
the following picture. In the quiescent phases the values 
of p s and R s fluctuate around constant values, yet their 
product remains invariant along the evolution of a halo. 

The S s invariance impli es that M s oc R 2 . This is very 
close to the M s oc R]. 92 of lZhao et all (120031) and the mi- 
nor discrepancy results from lZhao et al.l "not separating 
explicitly between the quiescent and violent phase s. We, 
therefore, support and validate one of I Zhao et al.l 's main 
points, namely the M s oc R" sc aling, and set a = 2. 
Thus, we validate also lZhao et aLl 's relation between the 
evolution of the concentration parameter and the MAH. 
Moreover, this relation can be easily extended to the 

5 II 

MAH — Q s relation, given the Q s Rs invariance. 

A final note concerns the entropy of a system of self- 
gravitating collisionlcss particles. The definition of the 
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Fig. 1. — The cosmological evolution of the main halo of the six models is shown as a function of the expansion parameter a. The violent 
phases are characterized by sudden increase (decrease) of R s , indicated by continuous lines, (Q s , dashed lines) and theses quantities remain 
roughly constant in the quiescent phase. At the violent phases the NFW fitting fails and the resulting R s and Q B become unstable. The 
thick horizontal lines correspond to the mean value of R s within the different quiescent phases. The WMAP3 model exhibits an extended 
initial phase of a succession of violent mergers that cannot be resolved into individual ones separated by quiescent phasess, like in the 
OCDM models. This is followed by a quiescent phase. The Q s curves have been shifted to fit within the window. 
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Fig. 2. — The cosmological evolution of the dimensionless phase-space density profile, Q(r/R s ), of the OCDMc (left panel) and the 
ACDM WMAP3a (right panel) halos. The Q(r/R s ) profile is plotted for a sample of snapshots for the two halos. The value of for each 
snapshot is indicated. The bottom panels show the the fractional residual deviation from a power law. The relative deviation from a power 
law behavior is smaller than twenty percent. The universal form of Q(r/R s ) is clearly exhibited. The Q(r/R s ) profiles of the other OCDM 
models arc virtually identical to the ones shown here. 



entropy (per particle) of a monoatomic ideal gas is givcin 

by 

s = k B ln(<3 -1 ) + const, (9) 

where ks is t he Bol tzmann constant (e.g ., 

iDalcanton fc Hoganl 12001k IWhite fc Naravanl Il987f ). 
Applying this definition to the DM particles provides 
one with a formal entropy of the DM. Our findings 
concerning structure and evolution of the phase-space 
density can be easily translated to the language of the 
entropy of the DM. It should be noted here that Eq. 



03 provides a local measure of the entropy. The long 
range nature of the gravitational interactions prohibits 
a simple extension of the entropy to an extensive 
quantity that characterizes the whole halo. Within this 
framework we can only refer to the entropy as a local 
property. 

6. DISCUSSION AND CONCLUSIONS 

The main elements of the structure and evolution of the 
DM halos can be summarized as follows. Halos evolve 
through two phases, quiescent and violent ones, which 
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represent the two extreme cases of smooth accretion and 
major mergers. In the quiescent phase the halo density 
distribution is closely approximated by the NFW pro- 
file. The inner (within R s ) halo mass surface density, 
S s , remains approximately constant throughout its en- 
tire evolution. Most importantly, the major mergers that 
take the halo from one quiescent phase to the other pre- 
serve the value of E s . During the quiescent phases the 
halo phase-space density profile follows a power law of 
the form Q(r) = Q a (r/R a ] r 13 with (3 « 1.9. The cosmo- 
logical evolution of the phase space density is given by 
Q s . Under the invariance of S s and the assumption of 
the virial equilibrium within R s , the evolution of Q s is 

dictated by R s , so that Q s oc R s ■ In the quiescent 
phases, Q s remains constant and it decreases discontin- 
uously in the violent phases. 

The Ti s invariance and the virial theorem (Eq. [6|) imply 
that the evolution of a halo proceeds while conserving its 
surface density and its energy per unit volume, or equiv- 
alcntly the pressure, within R a . The interesting point 
is that a typical halo undergoes a few violent events of 
major merging that destroy its equilibrium. Following 
each event it regains the NFW structure with a larger 
R Sl while preserving its pressure and the mean surface 
density within the new R s . Summarizing, a DM halo 
evolves along a sequence of NFW profiles, with an ever 
increasing R V i r and R s that grows only discontinuously, 



in the manner described in Paper I, while conserving the 
mean surface density and the pressure within R s . 

The evolution of the Q(r) profile of a given DM halo 
is predicted by Eq. [JJ The controlled numerical experi- 
ments of Papers I and II suggest that, for a given halo, 
T, s remains constant over a few quiescent phases inter- 
mitted by violent events. All the halos analyzed here 
(apart from WMAP3a) emerge from the same realiza- 
tion of the initial conditions that has been subjected to 
different constraints. In many ways these halos can be 
considered as a single halo that has been manipulated so 
as to modify its merging history. As such they all have 
roughly the same value of E s . The WMAP3a model pro- 
vides an independent realization of a DM halo in a differ- 
ent cosmology. The value of T, s is roughly twice as large 
than in all the other model. Of course the extremely 
poor statistics of our models cannot teach us about the 
scatter in T, s . The question arises as to what happens in 
the general case of a large ensemble of DM halos drawn 
from a large cosmological simulation. The scatter in E s 
is expected to determine the evolutionary tracks of the 
phasc-spasc density profiles. It is interesting to study 
the possible environment — T, s correlation, and to what 
extent this affects the evolution of the DM halos. 

The p s R s invariance implies that the lZhao et~afl (|2003f ) 
M s oc R™ scaling holds with a = 2, and so does the MAH 
- concentration parameter relation. It follows that this 
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Fig. 4. — The product E s = p s R s is plotted against the expansion parameter a for the six halo models. The panels present the individual 
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can be easily extended to calculate a similar MAH — Q s 



relation. The interesting ramification of this result is 



Evolution of phase-space density in dark matter halos 



7 



that any statistical algorithm for generating DM halos 
merging trees and/or MAHs can be extended so as to 
provide the evolution of the full NFW parameters of a 
halo along its merging tree. This can be used in semi- 
analytical modeling of galaxy formation in which realiza- 
tions of merging trees include the full NFW structure of 
the halos. 

The definition of Q(r) as a phase-space density is only 
a poor man's substitute for the 'real' phase-space den- 
sity in the six-dimensional phase space. Yet, pushing the 
analogy further on, a formal entropy is defined by Eq. [9] 
as a local variable. This local definition cannot be ex- 
tended to provide a global entropy of the halo. It follows 
that statements regarding the evolution of the entropy of 
DM halos are ill-defined, at least in the current context. 
Nevertheless, in the NFW scaling framework, in which 
the entropy of DM halos is taken as Q s , it increases with 
R s . In this context, the quiescent phases with R s ^const. 
correspond to adiabatic processes which preserve the en- 
tropy, and the violent phases — to non-adiabatic ones in 
which the entropy grows. 

The method of introducing an entropy by means of Eq. 
[5] can be considered as only formal, and its relation to a 
standard thermodynamic s needs to be questioned. Yet, 
Faltcnbac her et ail (|2007f ) have recently shown that the 
formally defined entropy of the DM is very closely as- 
sociated with the classical (ideal gas ) entropy of the in- 
tergalactic gas in clusters of galaxies. iFaltenbacher et all 
studied the entropy of the gas and DM in high reso- 
lution adiabatic SPH cosmological simulations — 'adia- 
batic' here is used in the sense that the entropy of the 
gas can only grow due to shock waves. Consequently, 
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following the accepted notation, the entropy of the gas 
can be expressed by means of K gas = ksT / p 2 / 3, , and 

by analogy Kdm = ^dm/Pdm ^ Q~ 2 ^ 3 - The power law 
behavior of Q{r) implies that Kdm cx r . Interestingly, 
outside the entropy core of r f=a i? s /2 the gas entropy fol- 
lows the DM, namely K gas oc Kom, and the calculated 
slope agrees well with observations. Faltcnbacher et al. 
amended the definition of K gas , such that the thermal en- 
ergy is extended to include the kinetic energy of the small 
scale (turbulent motion) and the DM and gas densities 
are normalized by their mean cosmological values. The 
newly defined K gas then coincides with the DM entropy 
outside the entropy core to a very good approximation. 
The analysis of Faltcnbacher et al. shows that, at least 
in galaxy clusters, the DM entropy coincides with the 
classical, ideal gas, entropy of the intracluster medium. 
It is not clear a priori why the gas and DM entropies 
should coincide, but numerical simulations suggest that 
they do. The simulations suggest that observational de- 
termination of the entropy of the intracluster gas should 
shed light and constrain the behavior of the DM entropy. 
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